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We present a quantum Monte Carlo study of the quantum correlations in the parametric lumines- 
cence from semiconductor microcavities in the strong exciton-photon coupling regime. As already 
demonstrated in recent experiments, a ring-shaped emission is obtained by applying two identical 
pump beams with opposite in-plane wavevectors, providing symmetrical signal and idler beams with 
opposite in-plane wavevectors on the ring. We study the squeezing of the signal-idler difference noise 
across the parametric instability threshold, accounting for the radiative and non-radiative losses, 
multiple scattering and static disorder. We compare the results of the complete multimode Monte 
Carlo simulations with a simplified linearized quantum Langevin analytical model. 
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In the last years, semiconductor microcavities in the 
strong exciton-photon couphng regimei"^ have been at- 
tracting a considerable deal of interest because of their 
remarkable nonlinear parametric interactions^'^i^i^'i: tak- 
ing advantage of a triply-resonant condition, ultra-low 
parametric oscillation thresholds have been observed 
in geometries which look very promising in view of 
applications. Very recently, experimental and theoret- 
ical investigations are starting to address the genuine 
quantum optical properties of the polariton parametric 
emissio n^i^d°d!d^d3,i4 ^ r^j^^ signal-idler pairs generated 
by the coherent scattering of two pump polaritons are 
expected to have non-classical properties, such as entan- 
glement and two-mode squeezing, which are interesting 
e.g. for quantum teleportation. The main limitation of 
the original non-degenerate parametric scheme where the 
cavity was pumped by a single incident beam at a finite 
'magic' angleSiiiii^, was the strong intensity asymmetry 
between the signal and idler photon emission. This signal- 
idler asymmetry is in fact strongly detrimental in view of 
the observation of significant extra-cavity quantum cor- 
relations to be used for continuous variable experiments. 

This difficulty has been overcome in recent 
experiments'^^ by using a pair of identical pump 
beams with small and opposite in-plane wavevectors. In 
this degenerate parametric scheme, a pair of perfectly 
symmetric signal and idler beams are emitted at the same 
frequency and with opposite wavevectors. For symmetry 
reasons, the momentum-space parametric luminescence 
pattern is in this case a ring, with approximately the 
same radius as the pump wavevector. Interestingly, this 
kind of ring-shaped polariton parametric luminescence 
can be obtained also with a single pump at normal 
incidence (zero in-plane wavevector) on a multiple mi- 
crocavity with multiple photonic branchesii. In order to 
quantify the performances of this system as a source of 
correlated photons, it is then important to characterize 
the robustness of the quantum correlations in the para- 
metric luminescence against competing effects such as 



radiative and non-radiative losses as well as multimode 
competition and multiple scattering processes. Given the 
unavoidable imperfections of any solid-state system, it is 
also crucial to assess the impact of a weak static disorder 
on signal-idler correlations: disorder is in fact known 
to be responsible for the so-called resonant Rayleigh 
scattering of pump photonsi^, which also produce a 
ring-shaped pattern in momentum space, yet without 
any quantum correlation. 

In this paper, we make use of the Wigner Quantum 
Monte Carlo methodic for polaritons in semiconductor 
microcavities to numerically tackle these key issues. The 
paper is structured as follows. In Sec. HI we present the 
model Hamiltonian and quantum Monte Carlo technique 
used to calculate the observables. Results for the ring- 
shaped polariton parametric luminescence with or with- 
out a static disorder are reported in Sec. HIl Correspond- 
ing numerical results for the quantum correlations are 
presented in Sec lIIII and then compared to a simplified 
quantum Langevin analytical model in Sec. IIVI Finally, 
conclusions are drawn in SecfVl 



I. HAMILTONIAN AND QUANTUM MONTE 
CARLO TECHNIQUE 

In this paper, we consider the quantum field Hamilto- 
nian introduced in Ref. S 

H = I ^ ^l(x)[h^^.-fy,(x)%]«'t(x) 

i3 = {X,C} 

+ dx*t^(x)*t^(x)#x(x)^'x(x) 

+ J dy:hFp{x,t)%{x) + h.c. , (1) 



where x is the in-plane spatial position. The field oper- 
ators "^x.ci^) respectively describe excitons and cavity 
photons. We assume an exciton density far below the 



2 



saturation density Ugair^, so the field operators obey the 



Bose commutation rules : [^'^(x), \I>J(x')] — 5^ ^i5ij. The 



linear Hamiltonian h° is 

'■J 



uJc{-iV) / ' 



(2) 



where wc(k) = w^-i/l + k^/A:^ is the cavity dispersion 
as a function of the in-plane wavevector k and kz is 
the quantized photon wavevector in the growth direc- 
tion. The exciton dispersion is assumed to be momentum- 
independent, i.e., wx(k) = ijj'^x- The quantity is the 
vacuum Rabi frequency of the exciton-cavity photon cou- 
pling. The eigenmodes of the linear Hamiltonian are 
called Lower and Upper Polaritons (iP, UP) . Their en- 
ergies are respectively hujLp{k) and fiwupik). The non- 
linear interaction term g is due to the exciton-exciton 
collisional interactions, which are modelled by a contact 
potentiali^. For the sake of simplicity, we restrict our- 
selves to the case of a circularly polarized pump beam, 
which allows us to ignore the spin degrees of freedom and 
the complex spin dynamics^SiSii^. The potential due to 
the static disorder is included in Vx,c(x). 

The polariton dynamics is studied by means of nu- 
merical simulations based on the so-called Wigner quan- 
tum Monte Carlo method, explained in detail in ReflTsl. 
Within this framework, the time-evolution of the quan- 
tum fields is described by stochastic equations for the C- 
number fields ^x.,c{^)- The evolution equation includes 
a non-linear term due to interactions, as well as dissi- 
pation and noise terms due to the coupling to the loss 
channels. Actual calculations are performed on a finite 
two-dimensional spatial grid of Ux x Uy points regularly 
spaced over the integration box of size x Ly. The dif- 
ferent Monte Carlo configurations are obtained as statis- 
tically independent realizations of the noise terms. 

Expectation values for the observables are then ob- 
tained by taking the configuration average of the stochas- 
tic fields. As usual in Wigner approaches'^, the stochas- 
tic average over noise provide expectation values for the 
totally symmetrized operators, namely: 



{Oi...On)w 



)...Op(jv)) , 



(3) 



the sum being made over all the permutations P of an 
ensemble of N objects. Each operator Oa represents 
here some quantum field component, while Oa is the 
corresponding C-number stochastic field. 

The relation between real- and momentum-space op- 
erators is: 



^'c(x) = ^i=^e^'^''ak, 



(4) 
(5) 



where fik (&k) represents the photonic (excitonic) de- 
struction operator for the k-mode, and satisfy the usual 



Bose commutation rules [ak,ajj.. 



[fek,^L 



3k, k' 

The expectation value of the in-cavity photon popula- 
tion fik = ajj flk in the k-mode reads: 

1 
2 



|ak| 



(6) 



where the overlined quantities are stochastic configura- 
tion averages, and ak is the C-number stochastic field 
value corresponding to the operator cik Because of the 
weak, but still finite transmittivity of the cavity mirrors, 
all observables for the in-cavity field transfer— i^l into the 
corresponding ones for the extra-cavity luminescence at 
the same in-plane momentum k. 



II. RESULTS FOR THE RING-SHAPED 
LUMINESCENCE 

A. In the absence of disorder 
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FIG. 1: Solid lines: energy dispersion of the Lower and Up- 
per Polariton branches. Dashed lines: the bare quantum well 
exciton and cavity photon dispersions. The filled circles in- 
dicate the wavevectors of the two pumps (pi and P2). Note 
that in the chosen configuration, the photonic fraction of the 
Lower Polariton at the pump wavevector is ~ 0.46. Cavity pa- 
rameters: WLr = 2.5 meV, hJ^c = 1400 meV, hw\ = 1400.5 
meV, fcz = 20/im~^. Pump parameters: kp — 0.6981 /im"^, 
ujj, — ajLp(kp) — 1398.2 meV. 



In this work, we will consider the following excitation 
field: 



(7) 



This field describes two identical monochromatic plane- 
wave pumps with opposite wavevectors oriented along 
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the X-axis. Both beams have the same values for the am- 
phtude Fp and the frequency ujp. This latter is chosen 
to be resonant with the LP-branch, i.e. LUp = a;ip(kp). 
Fig. [1] depicts the dispersion of the polariton branches 
and the position of the pump wavevectors. The scatter- 
ing process between a pair of ±kp pump polaritons via 
the non-linear interactions, gives rise to a pair of sig- 
nal/idler polaritons of opposite wavevectors ±k. Mod- 
ulo the weak blue-shift of the modes due to interactions, 
the energy-momentum conservation (phase-matching) is 
trivially fulfilled if Us — LUi = ojp and |k| = 
the |k| — kp parametric luminescence ring 
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FIG. 2: QMC results for the in-cavity photon population n^. 
Pump amplitude Fp/') = 5/^m~^ (just below the parametric 
instability threshold) . Number of Monte Carlo configurations 
: 330. The two rectangles denote the areas where the sig- 
nal and idler are integrated. Cavity parameters: hQ,R — 2.5 
meV, hio^ = 1400 meV, Hljx = 1400.5 meV, =20 /im"\ 
hflR = 2.5 meV, hjcx ^ hj ^ 0.1 meV, hg = 1.10"^ 
meV.^m"^. Pump parameters: kp — 0.6981 fj,m~^ , cop — 
uJLp(kp) ~ 1398.2 meV. Integration box size Lx = Ly = 90 
fim with rix ~ riy = 64 points; integration time step dt = 1.7 
fs. Using such a short time step has been necessary in order 
to obtain sufficient numerical precision on fourth-order field 
correlation functions. 



Fig. [5] shows the numerical results for the stationary 
state photon population inside the cavity for a value of 
the pump power below the parametric oscillation thresh- 
old: the ring-shaped parametric luminescence pattern is 
apparent. The interaction-induced blue shift of the po- 
lariton modes is responsible for the ring radius being 
slightly smaller than kp. Other interesting features can 
be observed in addition to the main ring: the strong 
spots at k = ±3kp are due to four-wave mixing processes 
(kp,kp) — > (±3kp,=Fkp); because of the stimulated na- 
ture of the underlying process, these spots fully inherit 
the coherence of the pump beams. Some luminescence is 
also observed along the x-axis in the vicinity of kp. Para- 
metric scattering processes involving polaritons from the 
same pump beam (kp,kp) — > (kp -I- 5k, kp — 5k) with 



|5k| ^ |kp| are responsible for this emission. As the pump 
beams are not tuned at the so-called magic angle, this 
emission is much weaker than the one on the ring. 

In the following, we will focus our attention on signal- 
idler pairs with wavevectors on the ring, and close to the 
y-axis {k^ — 0). To minimize discretization effects, we 
will average the signal/idler observables on the rectan- 
gular areas 2?s,i sketched in Fig. [21 which indeed contain 
quite a number of pixels. The corresponding photon pop- 
ulation operators Ns^i are defined as: 



(8) 



where Nd is the number of modes inside "Ds^ and fis^i 
are the average photon population operators. In term of 
the stochastic field, these latter read: 



(9) 



1- 

Q. 









i 




! - + - no disorder 




' —disordered potential 


, i'' 

r' 




r 





m ) 



FIG. 3: Left panel: QMC results for the density of the pump 
excitons pp = (bl^^bi^p) / (LxLy) . Right panel: signal/idler pho- 
ton populations na,i in the presence (solid line) or absence 
(dashed line) of the disordered potential. Same cavity and 
integration parameters as in Fig|2] 



The density of pump excitons Pp — {bj_j^ i>±\^^) / [LxLy) 
and the signal/idler populations ng^i are shown as a func- 
tion of pump power in the left and right panels of Fig. [3] 
respectively. As previously discusse d^^'^^ , the pump den- 
sity Pp smoothly increases up to the threshold for para- 
metric oscillation. The sub-linear dependence on power 
stems from the optical limiting effect due to the blue- 
shift of the ±kp modes by the repulsive interactions^. 
Around the threshold at i^p/7 ~ 5.75 /im~^, pp shows a 
downward kink, while the signal/idler populations have 
a sudden increase. For the realistic parameters used here, 
note how the density of excitons at the instability thresh- 
old remains moderate, and much lower than the exciton 
saturation density, pkp < 10^ cm~^ ^ psat- This shows 
the efficiency of the considered parametric process. 

In Fig. m we can see the temporal spontaneous build- 
up of the signal and idler luminescence starting from the 
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FIG. 4: QMC results for (hs) (Solid line) and (hi) (cir- 
cles) as a function of time (ps) for various pump amplitudes 
across the parametric oscillation threshold: -Fp/7 — 5 fim~^ 
(330 configurations), -Fp/7 = 5.5 nm~^ (340 configurations), 
Fp/j — 6/im^^ (180 configurations) and Fp/'y = 9/im~^ (60 
configurations). Same cavity and integration parameters as in 
FigEl 



vacuum fluctuations. For = the popula- 

tion of the signal/idler modes is still small Us^i ^ 1, 
while stimulated parametric scattering starts to be effec- 
tive for Fp/j = 5.5 /zm~^ when the occupation number is 
comparable or larger than 1. The parametric oscillation 
threshold has already been crossed for Fp/^ = 6/im~^. 
While the emission ring below threshold has the almost 
homogeneous intensity profile shown in Fig. [2l a symme- 
try breaking takes place above the threshold: a few modes 
are selected by mode competition effects, and a macro- 
scopic population concentrates into them, as shown in 
FigO It is interesting to note that that the number of 
Monte Carlo configurations needed for obtaining a given 
precision in the configuration average strongly depends 
on the regime under examination: as expected, much less 
simulations are required above the threshold. 



B. In the presence of static disorder 

The results in Fig. [21 [H and [5] have been obtained in 
the absence of static disorder, i.e. for Vc = Vx = 0. An 
arbitrary potential can be easily included in our calcu- 
lations. As a specific example, we have considered the 
disordered photonic potential reported in Fig. [SI consist- 
ing of a random ensemble of photonic point defects^^ . 
The corresponding emission pattern is shown in Fig. [7| 
for the same pump parameters as in the clean system of 
Fig. [21 The main effect of the disorder appears to be a 
significantly enhanced intensity on the luminescence ring. 
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FIG. 5: QMC results for the in-cavity photon population rik 
for a pump amplitude above the parametric oscillation thresh- 
old Fp/'y = 6 ^.Ta"^ . Number of Monte Carlo configurations : 
180. Same cavity and integration parameters as in Figl2] 
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FIG. 6: Disordered photonic potential (in meV/ft) used for 
the simulations of Fig. I3I7I and (TO] 



This occurs because of the resonant Rayleigh scattering 
of each of the pumps. Note also the weak "eight" -shaped 
pattern— due to the parametric amplification of the res- 
onant Rayleigh scattering ring. 

The signal population as a function of the pump power 
is plotted in the right panel of Fig. [31 below threshold, 
the photon population in the presence of disorder is much 
larger than in the clean system. On the other hand, the 
difference between the two populations is much less im- 
portant above threshold when the non-linear stimulated 
parametric scattering dominates over the linear Rayleigh 
scattering processes. Despite the very different low inten- 
sity behaviour, the threshold is reached in both cases at 
values close to Fp/7 ~ 5.5 ^m~^. 
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III. QUANTUM CORRELATIONS 



FIG. 7: QMC result for the in-cavity photon population rik 
in the presence of the disordered potential shown in Fig. [U] 
The image saturates the gray scale. Pump amplitude Fp/'y — 
5/im^^. Number of Monte Carlo configurations : 50. Same 
cavity and integration parameters as in Fig[2] 



In the present section we study the correlation prop- 
erties of the signal and idler emissions. For the sake of 
simplicity, we restrict our attention here to those fluc- 
tuations which are associated to the intrinsic losses of 
the parametrically emitting system, and we neglect all 
other possible noise sources that may appear in actual 
experimental setups, e.g. pump intensity fluctuations. To 
characterize the quantum nature of the correlations^^ be- 
tween the signal and idler modes, it is useful to consider 
the quantity N± = Ns±Ni i.e. the sum and difference of 
the signal and idler photon population. The correspond- 
ing normalized noise reads: 



(Nl) ~ (iV±)2 _ (^s) - (^.>' + (Nf) - (iV,)2 ± 2 ((iV.iV,) - (TV,) (TV,)) 



(10) 



Hence, the fourth-order moments of the fields (NjNk) — X^keD ic'eCk (""^^k') (where j,k £ {s,i}) play a key-role in 
the determination of the quantum behaviour of the system. In terms of the averaged stochastic quantities we have: 



(iVj) 



- ^ /akflkakak + ajj.akajj.ak + ajj.akakajj. -I- akajj.ajj.ak + akajj.akdjj. -1- akakajj.ajj.\ - ^(aj(.ak -|- Okajj.) 



kei?,- 



{nk){nk') - (^k)(fik) 




E 



1 



(11) 



with j = (i, s). The intensity correlation between signal and idler modes is: 



(NsNi) = ^ E (("k^k + akSjj. - 1) (ajj.,ak' + Sk'ajj., - i; 

keVs 
k'eVi 

= E E 



keVs k'eVi 

For uncorrelated and shot-noise limited signal and idler 
beams, one would have cr^ = 1: this value is the so- 
called Standard Noise Limited. Having cr~ < 1 means 
that non-classical correlations exist between signal and 
idler, in particular a squeezing of the difference intensity 
noise21i22i. As the polariton states are half-photon half- 



(12) 




k'eVi 



exciton, the optimal noise reduction of the photon field is 
reduced by half with respect to an ideal x^^^ purely pho- 
tonic parametric oscillator system-S; noise reduction does 
not concern in fact the photon field taken independently, 
but rather the whole polariton field. As long as we neglect 
multiple scattering and disorder effects, this is the main 
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difference compared to standard x*^^) parametric oscilla- 
tors; an analytical model for these issues will be provided 
in SeclIVI Note that throughout all the present paper 
we are interested in one-time correlations; the difference 
noise is therefore integrated over all the frequencies and 
no frequency filtering is considered. Some frequency fil- 



tering around uj 



would purify the squeezing as the 



quantum correlations are largest around this value of the 
frequency's. 
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FIG. 8: Left panel: QMC results for the time-evolution of the 
normalized photonic sum and difference noises a* in the ab- 
sence of disorder. Right panel: in-cavity photon population nk 
at t = 1800 ps. Pump amplitude -Fp/7 = 5 ^ra~^ . Number of 
Monte Carlo configurations: 330. Same cavity and integration 
parameters as in Fig[2] 



In Fig. [5] we have plotted the time dependance of the 
normalized noises cr^ for a given Monte Carlo realization 
and in the absence of disorder. The stationary-state av- 
erage values for the same quantities are plotted in Fig. 
ini as a function of the pump intensity. Quantum correla- 
tions cr^ < 1 exist in the difference noise at low intensi- 
ties, while it monotonically increases towards — 1 for 
higher intensities, making the signal/idler correlations al- 
most purely classical well above threshold^l. No specific 
feature is found in this quantity at the threshold. 

On the other hand, the sum noise is always above 
the standard noise limit, and shows a sudden increase 
at the parametric threshold. The fact that well above 
the threshold it does not go back to the standard noise 
is due to the presence of several competing parametric 
oscillation modes. Depending on whether the oscillating 
modes lay inside or outside the selected regions Vg^i, the 
signal/idler populations n^.i vary between and their 
maximum value, while remaining almost equal to each 
other. This implies that the sum noise is large, of 
the order of the signal/idler populations while the 
difference noise a~ remains small. 

In Fig llOi we have analyzed the sum and difference 
noise in the presence of a disordered potential. For the 
same value of pump intensity, the difference noise a~ 
is now somehow larger than in the absence of disorder: 
the resonant Rayleigh scattering creates in fact unpaired 
photons into the luminescence ring and deteriorates the 
pair correlations between the signal and the idler. For 
very low intensities, the dominant contribution comes 
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FIG. 9: QMC results for the normalized photonic sum and 
difference noises cr* as a function of the pump power in the 
absence of disorder. The difference noise is fitted by a poly- 
nomial curve, while the line in the plot for is a guide for 
the eye. Same cavity and integration parameters as in Fig[2] 



from the Rayleigh scattering processes implying that 
both the sum and the difference noises have to tend 
toward the Standard Noise Limit — \. Because of the 
competition between the Rayleigh and the parametric 
scattering, the difference noise a~ attains its minimum 
in the vicinity of the threshold and then increases 
because of the increasing noise of the two beams. As 
disorder is able to mix the modes respectively inside and 
outside the selected regions I^s^i, the difference noise 
can grow above a~ = 1 at high pump powers. For the 
same reason, the sum noise cr"*" has a weaker growth 
above the threshold than in the absence of disorder. This 
physical interpretation of the role of the disorder has 
been confirmed by several other simulations (not shown) 
performed with different realizations of the disordered 
potential. 
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FIG. 10: QMC results for cr^ as a function of the pump in- 
tensity in the presence of disorder. The lines are a guide for 
the eye. Same cavity and integration parameters as in Fig[2] 
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IV. SIMPLIFIED ANALYTICAL MODEL 

The aim of this section is to compare the resuhs of the 
complete Quantum Monte-Carlo calculations to a simpli- 
fied input-output analytical model based on a lineariza- 
tion of the Hamiltonian. This is done by treating the 
intense pump as a classical, undepleted, field, i.e. replac- 
ing the pump mode operators with their mean-field ex- 
pectation values. Obviously, this approximation is valid 
only well below the parametric oscillation threshold. Con- 
centrating our attention on those processes which satisfy 
the phase-matching condition and neglecting all the non- 
resonant others, we can write the linearized Hamiltonian 
in the following, simplified form: 



hflji{ fej^flk -I- a^^k 



, (13) 



where 6k is the exciton creation operator, ojx = + 
L^L (I'^iP + l'P2p) is the blue-shifted exciton frequency 
because of interactions, and k — ^ 'Pi 7^2 is the ef- 
fective parametric interaction constant in terms of the 
pump fields 7'i_2e~''^p* = (6±k {t))- Taking the standard 
vacuum as the initial state of the photon and exciton 
fields, the expectation values of the quantum Langevin 
forces are: 



{FcavA^Wl^^M) = 4^r,,,M<5(cc;-c^')'5k,k' ,(14) 
{Fexc.M['^]Fl^ck'i'^]) = 47rrea:cM'5(a; - w')<5k,k' ,(15) 



where Tj[Lu] is the complex broadening due to the cou- 
pling to the external bath. Since the relevant spectral 
domain in the degenerate parametric process is con- 
centrated around cjp, we are allowed to simplify the 
treatment by taking frequency independent linewidth 
^ph,exc['^] = 7c,x/2. The quantum Langevin equations 
in frequency space read: 



I OkM 

6k [w] 



with the matrix Alk.cj.cj defined for i = X^C as: 



-Pcai;,k[w] ^ 
-Fea;c,k[^] 

^cL,-k[2^p-^] 

V^eL-k[2^p-H } 



(16) 



My 



-Ac{uj-2ojp)-i-fc/2 

-K* 





K 



in terms of Aj(a;) — cui ~ to . 

The relation between the time dependent and fre- 
quency dependent operators is: 



which leads to: 



-Ax{uj~2ujp)^i-fx/2/ 



ak(<) = / ^flkMe" 

J ^TT 



(18) 



(iV,)= ^ / |^(7c|ei3rH+7x|ei4|'H) 



where ak[w] is the component at ui of the photonic de- 
struction operator for the k-mode. In the following, we 
will call Gk,ui,ujp = ^*-^k;Li.jp- 1"^^ signal photon pop- 
ulation operator ns{t) inside the cavity can be written 
as: 



(17) 



(20) 



duj duj' 



(19) 



To calculate the sum and difference noise, the second 
order momenta (7y2) _ (^Nsf, (Nf) - (Ni)^ and (N,Ns) + 
(NsNi) — 2{Ni){Ns) are needed. After some algebra, we 
obtain the final expressions (for j = {s,i)): 



duj' 



i3p+7xiei4p)H / ^(7c|t;iiP+7xiei2p)K] 



(21) 



dLUlduj2 

'WW 



E ('^C'SiJk, 2ujp - uji]gi^[-k,uji] + 7x^12 [k,2wp - uji]gi^[~k,uJi]) 



(7c^/ii[k, 2wp - W2]5i3[-k,cj2] +7x012 [k, 2tjp - Lj2]^i4[-k,W2]) . (22) 
The results for and cr+ are plotted in Fig. [TT] The qualitative similarities between these results and the 
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cise estimation than the QMC calculation. As we have 
already mentioned, partition noise due to the half pho- 
ton half exciton nature of the polaritons is responsible for 
a significantly larger value of cr" than in standard x^^-* 
parametric emitters^^. While in the QMC calculations 
the inclusion of static disorder has been done straight- 
forwardly, the simplified analytical model can not be ex- 
tended easily to the disorder case, and, most of all, it 
would lose all its simplicity. 



Pump exciton density (n m 



Pump exciton density (n m 



FIG. 11: Analytical results for the photonic difference and 
sum noises cr* as a function of the pump exciton density pp 
in the absence of disorder. Same cavity and integration pa- 
rameters as in Fig[2] hujp — hujLp(k.p) + 0.07 meV = 1398 
meV and huj°x = 1400.1 meV. 



QMC ones (without disorder) of Fig [5] are apparent for 
pump intensities up to the parametric threshold. Here, 
the linearized model breaks down, as it predicts a di- 
verging signal/idler intensity. Although the predictions 
for the threshold pump intensity differs from the QMC 
one by approximately twenty percent, still the analytic 
value in the low intensity limit a~ ~ 0.75 is well within 
the (quite large) error bars of the QMC simulations with- 
out disorder—. In this limit, the analytical calculation, 
that neglects interactions between the signal and idler 
modes, becomes indeed exact and provides a more pre- 



V. CONCLUSION 

In conclusion, wc have presented a quantum Monte 
Carlo study of the quantum correlations in the ring- 
shaped parametric luminescence from semiconductor mi- 
crocavities in the strong exciton-photon coupling regime. 
Our results suggest that even in presence of multiple scat- 
tering, realistic losses and static disorder, the signal and 
idler beams maintain a significant amount of quantum 
correlations. The dependance of the intensity quantum 
correlation on the pump intensity has been characterized 
across the parametric instability threshold, showing the 
regime where the non-classical features are maximized. 

We thank G. Bastard, J. Bloch, A. Bramati, C. 
Diederichs, E. Giacobino, J-Ph. Karr, C. Leyder, N. Reg- 
nault, M. Romanelli, Ph. Roussignol and J. Tignon for 
discussions. 



* E-mail: 'cristiano.ciuti@univ-paris-diderot.fr| 

^ J. Baumberg and L. Vifia (Eds.) Special issue on Micro- 
cavities, [Semicond. Sci. Technol. 18, S279 (2003)]. 

^ B. Deveaud (Ed.), Special issue on the "Physics of semicon- 
ductors microcavities" , Phys. Stat. Sol. B 242, 2145-2356 
(2005) and references therein. 

^ J. J. Baumberg, P. G. Savvidis, R. M. Stevenson, A. I. 
Tartakovskii, M. S. Skolnick, D. M. Whittaker and J. S. 
Roberts, Phys. Rev. B 62, R16247 (2000) 

* C. Diederichs, J. Tignon, G. Dasbach, C. Ciuti, A. 
Lemaitre, J. Bloch, Ph. Roussignol and C. Delalande, Na- 
ture 440, 904 (2006). 

^ P. G. Savvidis, J. J. Baumberg, R. M. Stevenson, M. S. 

Skolnick, D. M. Whittaker and J. S. Roberts, Phys. Rev. 

Lett. 84, 1547 (2000). 
® G. Dasbach, M. Schwab, M. Bayer, D.N. Krizhanovskii and 

A. Forchel, Phys. Rev. B 66, 201201 (R) (2002). 

G. Dasbach, A. A. Dremin, M. Bayer, V. D. Kulakovskii, 

N. A. Gippius and A. Forchel, Phys. Rev. B 65, 245316 

(2002). 

* I. Carusotto and C. Ciuti, Phys. Rev. Lett. 93, 166401 
(2004). 

® J. Ph. Karr, A. Baas, R. Houdre and E. Giacobino, Phys. 

Rev. A 69, 031802(R) (2004). 
^° 3. Ph. Karr, A. Baas and E. Giacobino, Phys. Rev. A 69, 

063807 (2004). 
" C. Ciuti, Phys. Rev. B 69, 245304 (2004). 



S. Savasta, O. Di Stefano, V. Savona and W. Langbein, 
Phys. Rev. Lett. 94, 246401 (2005). 

I. Carusotto and C. Ciuti, Phys. Rev. B 72, 125335 (2005) 

A. Baas, J.-Ph. Karr, M. Romanelh, A. Bramati and E. 

Giacobino, Phys. Rev. Lett. 96, 176401 (2006). 

P. G. Savvidis, C. Ciuti, J. J. Baumberg, D. M. Whittaker, 

M. S. Skolnick and J. S. Roberts, Phys. Rev. B 64 075311 

(2001). 

^® M. Romanelli , C. Leyder, J-Ph. Karr, E. Giacobino and 
A. Bramati, Phys. Rev. Lett. 98, 106401 (2007). 
C. Diederichs, C. Leyder, A. Bramati, D. Taj, T. Lecomte, 
A. Lemaitre, L. Largeau, J. Bloch, C. Ciuti, Ph. Roussig- 
nol, C. Delalande, E. Giacobino, J. Tignon, unpublished. 
See, e.g., H. Stolz, D. Schwarze, W. von der Osten, G. 
Weimann, Phys. Rev. B 47, 9669 (1993); R. Houdre, C. 
Weisbuch, R. P. Stanley, U. Oesterle, M. Ilegems, Phys. 
Rev. B 61, R13333 (2000); W. Langbein and J. M. Hvam, 
Phys. Rev. Lett. 88, 047401 (2002). 

C. Ciuti, P. Schwendimann and A. Quattropani, Semicond. 
Sci. Technol. 18, S279 (2003). 

A. Kavokin, G. Malpuech and M. Glazov, Phys. Rev. Lett. 
95, 136601 (2005). 

21 I. Shelykh, G. Malpuech, K. V. Kavokin, A. V. Kavokin 
and P. Bigenwald, Phys. Rev. B 70, 115301 (2004). 

22 P. G. Lagoudakis, P. G. Savvidis, J. J. Baumberg, D. 
M. Whittaker, P. R. Eastham, M. S. Skolnick and J. S. 
Roberts, Phys. Rev. B 65, 161310(R) (2002). 



9 



L. Mandel and E. Wolf, Optical Coherence and Quantum 
Optics (Cambridge University Press, Cambridge, 1995). 
D.F. Walls and G.J. Milburn, Quantum Optics (Springer- 
verlag, Berlin, 1994). 

M. Wouters and I. Carusotto, Phys. Rev. B 75, 075332 
(2007). 

We have verified that the results do not qualitatively 
change when difi'erent specific potential configurations are 
used. 

N. Treps and C. Fabre, L aser Phys. 15, 187 (2005), 
( Preprint quant-ph/0407214') . 

C. Fabre, S. Reynaud in Quantum fiuctuations. Lecture 
Notes of Les Houches Summer School Session LXIII, edited 
by S. Reynaud, E. Ciacobino and J. Zinn- Justin. (North- 
Holland 1995). 

J. Laurat, T. Coudreau, N. Treps, A. Maitre, and C. Fabre, 
Phys. Rev. Lett. 91, 213601 (2003). 

C. Fabre in Confined Photon systems. Fundamentals and 



Applications, Lectures Notes in Physics of Cargese Sum- 
merschool, edited by H. Benisty, J-M Gerard, R. Houdre, 
J. Rarity and C. Weisbuch (Springer 1998). 
^1 C. Ciuti and I. Carusotto, Phys. Stat. Sol. B 242, 2224 
(2005). 

C. Ciuti, P. Schwendimann, B. Deveaud and A. Quat- 
tropani, Phys. Rev. B 62, R4825 (2000). 

M. Saba, C. Ciuti, J. Bloch, V. Thierry-Mieg, R. Andre, 
Le Si Dang, S. Kundermann, A. Mura, G. Bongiovanni, 
J. L. Staehli and B. Deveaud, Nature (London) 414, 731 
(2001). 

D. M. Whittaker, Phys. Rev. B 63, 193305 (2001). 

C. Ciuti and L Carusotto, Phys. Rev. A 74, 033811 (2006). 
^® A straighforward application of squ eezing theory shows 
that the exact value a~ = 0.75 is recovered in the limit 
where the upper polariton can be neglected and the lower 
polariton is exactly half photon and half exciton. 



